A finite element strategy coupling a gradient-enhanced damage model and cohesive cracks for quasi-brittle materials

A new combined strategy to describe failure of quasi-brittle materials is presented thus allowing the complete description of the process, from initiation of damage to crack propagation. For the early stages of the process, and in order to overcome the well-known problems characterising local descriptions of damage (e.g. mesh-dependence), a gradient-enhanced model based on smoothed displacements is employed. In order to deal with material separation, this continuous description is coupled to a cohesive crack when damage parameter exceeds a critical value. Some difficulties may arise when dealing with the transition from regularised damage models to evolving cracks: crack initiation, crack-path direction, energetic equivalence... In this work, a discrete cohesive crack is introduced when the damage parameter exceeds a critical value. On the one hand, and to determine the crack-path direction, the medial axis of the already damaged profile is computed. That is, a geometric tool widely used in the computer graphics field is used here to track the crack surface. Since this technique is exclusively based on the shape of the regularised damage profile, no mesh sensitivity is observed when determining the crack direction. On the other hand, and to define the cohesive law, an energy balance is imposed thus ensuring that the fracture energy not yet dissipated in the damage zone is transferred to the crack

Preface

Peer ReviewedPostprint (published version

Accuracy of two stress update algorithms for shear-free large deformations paths

The behavior of two stress update algorithms for shear-free large deformation paths is analyzed. The first algorithm has a truncation error of order 1. The second algorithm has a truncation error of order 2. As a consequence, the global performance of the second algorithm is clearly superior. However, for the particular case of shear-free deformation paths, the first algorithm correctly predicts null shear stresses, while the second one does not. This behavior was reported in a previous paper for an extension-rotation test. In this note a general shear-free deformation path is considered in full detail.Peer ReviewedPostprint (author’s final draft

Accuracy of two stress update algorithms for shear-free large deformations paths

The behavior of two stress update algorithms for shear-free large deformation paths is analyzed. The first algorithm has a truncation error of order 1. The second algorithm has a truncation error of order 2. As a consequence, the global performance of the second algorithm is clearly superior. However, for the particular case of shear-free deformation paths, the first algorithm correctly predicts null shear stresses, while the second one does not. This behavior was reported in a previous paper for an extension-rotation test. In this note a general shear-free deformation path is considered in full detail

A Comparison of two objective stress rates in object-oriented codes

An object-oriented code, CASTEM2000, has provided the framework for algorithmic development and numerical testing. In these codes, information is stored and manipulated as objects which represent the entitiesrelevant to a finite element method computation

Numerical analysis of nonlinear large-strain consolidation and filling

Finite strain consolidation and filling of soft sediments at high water level is a challenging problem because of its highly non-linear physical and mathematical aspects. Several numerical schemes designed for this problem are presented as well as simple numerical improvements for a better handling of the extremely high variations of the material properties with depth. The numerical algorithms developed are robust and verify convergence of the iterative schemes instead of the more classical approaches based on choosing time increments ‘sufficiently’ small and assuming convergence at every step. A set of computer programs has been developed to predict magnitude and rate of large-strain self-weight one-dimensional and pseudo bi-dimensional (i.e. one-dimensional deformation, bi-dimensional flux) consolidation during and after deposition, that is, coupling filling and consolidation phenomena. The actual life of the deposit can be numerically simulated combining filling periods and quiescent periods where surcharges (or capping) can exist. Consequently, they are a basic technique for the design of disposal ponds

Hygroscopicity issues in powder and grain technology

The effect of hygroscopicity on flowability of powders and bulk solids with applications in the packaging industry is experimentally and numerically investigated. Firstly, four granular materials are tested at different water contents to study the impact of relative humidity on some hydro-mechanical properties, namely the hydraulic diffusivity on wetting, as well as the shear strength and compressibility properties of the materials. Next, a capillary model covering a wider water content range –compared to the previous tests– is applied to discrete element simulations of a granular column collapse set-up. These simulations give further insight into important aspects of grain hygroscopicity in packaging and other industrial applications (such as the kinematics of flow and flowability issues), which are outside the scope of conventional experimental testing.Postprint (published version

Two stress update algorithms for large strains: accuracy analysis and numerical implementation

Two algorithms for the stress update (i.e., time integration of the constitutive equation) in large-strain solid mechanics are compared from an analytical point of view. The order of the truncation error associated to the numerical integration is deduced for each algorithm a priori, using standard numerical analysis. This accuracy analysis has been performed by means of a convected frame formalism, which also allows a unified derivation of both algorithms in spite of their inherent differences. Then the two algorithms are adapted from convected frames to a fixed Cartesian frame and implemented in a small-strain finite element code. The implementation is validated by means of a set of simple deformation paths (simple shear, extension, extension and compression, extension and rotation) and two benchmark tests in non-linear mechanics (the necking of a circular bar and a shell under ring loads). In these numerical tests, the observed order of convergence is in very good agreement with the theoretical order of convergence, thus corroborating the accuracy analysis

Efficient and reliable nonlocal damage models

We present an efficient and reliable approach for the numerical modelling of failure with nonlocal damage models. The two major numerical challenges––the strongly nonlinear, highly localized and parameter-dependent structural response of quasi-brittle materials, and the interaction between nonadjacent finite elements associated to nonlocality––are addressed in detail. Reliability of the numerical results is ensured by an h-adaptive strategy based on error estimation. We use a residual-type error estimator for nonlinear FE analysis based on local computations, which, at the same time, accounts for the nonlocality of the damage model. Efficiency is achieved by a proper combination of load-stepping control technique and iterative solver for the nonlinear equilibrium equations. A major issue is the computation of the consistent tangent matrix, which is nontrivial due to nonlocal interaction between Gauss points. With computational efficiency in mind, we also present a new nonlocal damage model based on the nonlocal average of displacements. For this new model, the consistent tangent matrix is considerably simpler to compute than for current models. The various ideas discussed in the paper are illustrated by means of three application examples: the uniaxial tension test, the three-point bending test and the single-edge notched beam test